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ABSTRACT 

We present a complete sample of 29 7-ray bursts (GRBs) for which it has been possible 
to determine temporal breaks (or limits) from their afterglow light curves. We interpret these 
breaks within the framework of the uniform conical jet model, incorporating realistic estimates of 
the ambient density and propagating error estimates on the measured quantities. In agreement 
with our previous analysis of a smaller sample, the derived jet opening angles of those 16 bursts 
with redshifts result in a narrow clustering of geometrically-corrected 7-ray energies about = 
1.33 X 10^^ h^'^ erg; the burst-to-burst variance about this value is 0.35 dex, a factor of 2.2. 
Despite this rather small scatter, we demonstrate in a series of GRB Hubble diagrams, that 
the current sample cannot place meaningful constraints upon the fundamental parameters of 
the Universe. Indeed for GRBs to ever be useful in cosmographic measurements we argue the 
necessity of two directions. First, GRB Hubble diagrams should be based upon fundamental 
physical quantities such as energy, rather than empirically-derived and physically ill-understood 
distance indicators (such as those based upon prompt burst time-profiles and spectra). Second, a 
more homogeneous set should be constructed by culling sub-classes from the larger sample. These 
sub-classes, though now first recognizable by deviant energies, ultimately must be identifiable by 
properties other than those directly related to energy. We identify a new sub-class of GRBs ( "f- 
GRBs") which appear both underluminous by factors of at least 10 and exhibit a rapid fading 
(/,y oc t^^) at early times (t <; 0.5 day). About 10-20% of observed long-duration bursts appear 
to be f-GRBs. 

Subject headings: cosmology: miscellaneous — cosmology: observations — gamma rays: bursts 
— distance scale — supernovae: individual (1998bw, 2003dh) — cosmological parameters 



1. Introduction 



The observational establishment of the extra-Galactic origin of gamma-ray bursts (GRBs) (Metzger 
et al. 1997) enlivened the possibility of exploiting GRBs for cosmological investigations. The discovery of 
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GRBs at high rcdshifts (GRB 971214, z = 3.42, Kulkami et al. 1998; GRB 000131, z = 4.51, Andersen 
et ah 2000) further demonstrated that GRBs, Uke quasars, could be used to study the early Universe and in 
particular the intergalactic medium, the interstellar medium of young galaxies and star-formation in dusty 
regions (see Loeb & Barkana 2001 for a review). 

The developments in the GRB field thus appear to mirror the early years of quasar astronomy. In both 
cases, astronomers were astonished by the unexpectedly large brilliance. In both cases, the fields split into 
two camps: many who used the objects for cosmology and some who were intrigued by the workings of the 
central engine. 

However, the large luminosity function of quasars precluded their use as a means to determine the 
geometry of the Universe, i.e. for cosmography. The bewildering diversity of GRBs in almost any respect - 
fluence, peak luminosity - likewise seemed to make them unsuitable for this purpose. 

Frail et al. (2001; hereafter FKS) presented a complete sample of GRBs with good afterglow follow-up 

and known rcdshifts. Most afterglows exhibit a steepening, seen either at optical or radio wavelengths, which 
is expected if the explosion geometry is conical (opening angle, 0jet) rather than spherical (Rhoads 1997a; Sari 
et al. 1999). FKS interpret the breaks making several simplifying assumptions (constant ambient density, 
constant 7-ray conversion efficiency) and find the surprising result that the seemingly most energetic bursts, 
E\so{^l) (isotropic-equivalent 7-ray energy release), possess the smallest beaming factor, /{, = ! — cos(0jet)- 
In that sample, the true energy release, = fbEisoi^), was strongly clustered. 

This finding - that GRBs are standard-energy explosions - opens up the possibility of using GRBs and 
their afterglows for cosmography. Indeed, in the August 2001 meeting at Garching^ we noted that a Hubble 
diagram constructed from in the FKS sample had a smaller scatter (0.31 dcx = 0.79 mag) compared to 
the dust-corrected _B-band peak magnitude of Type la supernovae (SNc) as of 1992 ((j=0.84mag; see Figure 
3 of van den Bergh & Pazder 1992). Recently, Schaefer (2003) presented another GRB Hubble diagram based 
on bursts that have two empirical distance indicators: the variability index (Fenimore & Ramirez-Ruiz 2003; 
Norris ct al. 2000; Rcichart ct al. 2001) and the lag- luminosity correlation (Norris 2002). The Schaefer 
analysis find a similarly small scatter about a standard candle. The current distance-indicator sample is 
small (9 sources) and restricted to GRBs with known redshifts and 7-ray observations from BATSE. 

As the experience with previous efforts at cosmography has shown, standard candles whose physics 
is understood (e.g., Cepheid variables) have proved more useful and robust than those which were based 
on poorly-understood or unphysical phenomena (e.g., brightest stars, number counts). For Type la SNe 
the physics of the explosion is reasonably well understood. Even so, for Type la SNe, it is still important 
to understand the intrinsic diversity (i.e., the distinction between the so-called sub-classification of "Branch 
normals" and "Branch peculiars" ; Branch et al. 1993), to eliminate the large fraction (-^36%; Li et al. 2001) of 
outliers on a sound empirical basis, and thereby increase the precision of the standard candle. Ultimately, the 
distance indicator must be calibrated to high enough accuracy to be competitive with precision cosmological 
experiments such as WMAP (Spcrgel et al. 2003). 

As with Type la SNe, the apparent constancy of the prompt energy release in GRBs is plausible on 
physical grounds. Furthermore, indications of a constant energy release come from studies of GRB afterglows 

(Piran et al. 2001; Panaitescu & Kumar 2002; Berger et al. 2003), a physically distinct emission process from 
th{^ prompt biu'st (unission. It is rc^asonablc to expect that complementary afterglow data (as with photometry 



^http: / /www.mpa-gaTching.mpg.de/~light/ 
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and spectroscopy of Type la SNe) would allow us to identify potential families, or sub-classes, in GRBs and 

thus compile more homogeneous samples. 

With this motivation and noting that the sample of GRBs with quality afterglow data has increased 
since the FKS analysis (from 17 to 29), we revisit the current state of the GRB Hubble diagram paying 
careful attention to the error analysis and identifying potential sub-classes of GRBs. 

2. Physical interpretations for observed temporal breaks 

Temporal breaks (scciilar steepening of a light curve) have been observed or inferred at optical, X-ray 
and radio wavelengths in a number of burst afterglows. When well-sampled broadband imaging exists, these 
breaks appear to be achromatic or very nearly achromatic. The (near) achromaticity is suggestive of a 
geometric rather than an emissivity effect (such as expected from an hydrodynamically evolving synchrotron 
shock). 

At present, two different geometric pictures are used to interpret the existence of ibrcak, the characteristic 
break time between two asymptotic power- law decays (see Kumar & Granot 2003 for a comprehensive review). 
In one interpretation (the "top hat" model), the energy per unit solid angle is constant within some cone 
with a half-angle opening Oj^t and falls to zero outside this angle. The break occurs at a time when the 
bulk Lorcntz factor of the shock has slowed to F « (Rhoads 1997b; Sari et al. 1999). In the other 

interpretation (the "universal jet" or "structured jet" model), the energy per unit solid angle varies across 
the face of the jet and the observer lies at an angle ©obs from the jet axis (Meszaros et al. 1998; Zhang & 
Meszaros 2002). Here, the break occurs at a time when the bulk Lorcntz factor of the shock in the direction 
of the observer slows to r(9obs) ~ l/r(0obs)- The present data cannot yet distinguish between the top-hat 
and universal jet models (Rossi et al. 2002). The observational implication of a hybrid of the two models 
has only recently been explored (Kumar & Granot 2003). 

The inference of 6jet or 6obs from tbreak requires an estimate of the evolution of F as a function of 
observer time. The dynamical evolution of F is theoretically prescribed by the relativistic Sedov solution 
(Blandford & McKec 1976). This solution is dependent upon the isotropic-equivalent energy in the shock 
and the nature of the ambient medium into which the shock propagates for t < tbreak (e.g.. Sari et al. 
1998). The isotropic-equivalent energy in the prompt burst is measured using the observed fluence, prompt- 
burst spectrum, and distance to the source. The nature of the ambient medium is generally assumed to be 
homogeneous — density as a function of distance r from the explosion site as p{r) = nor° — or wind-stratified 
{p{r) = A*r~^). Light curves of the more general case p{r) oc have also been explored (e.g., Panaitescu 
et al. 1998). 

Of the ~30 well-studied GRB afterglow to date, there has been only f«2 bursts for which a wind- 
stratified media are strongly favored statistically over homogeneous media (GRB 980519: Frail et al. 2000a 
and GRB 011121: Price et al. 2002b). In the remainder, winds are either excluded (e.g., 990510, Harrison 
et al. 1999; Panaitescu & Kumar 2002; 990123 Chevalier & Li 2000) or allowed at similar significances as 
homogeneous profiles (e.g., 980329; Yost et al. 2002; see also Panaitescu & Kumar 2002). Recently, Panaitescu 
& Kumar (2003) have reanalysed the afterglow of GRB 990510 and showed an improved consistency with 
the wind hypothesis by taking into account a structured jet. 

In what follows, we construct a set of geometry-corrected energies for GRBs, following from the formalism 
of FKS, under the assumption that the early afterglow propagates into a density of p{r) = nor° (we later 



-4- 



address the differences when, instead, a wind-stratified medium is assumed). We also assume a top-hat model 
for the structure of the jet but note the resulting energy determinations are identical under the universal jet 
model interpretation when the energy per unit solid angle drops as (Zhang & Meszaros 2002; see also 
Kumar & Granot 2003 and referenced therein). 



3. Formalism and Error Analysis 



Following the definitions and notation of FKS and Bloom et al. (2001b), the value for Ej, the total 
prompt energy release in a certain bandpass, may be found from: 



S i^fcf. 



(1) 



where S-y is the fluence received in some observed bandpass and Di is the luminosity distance at redshift z. 
The quantity A; is a multiplicative correction of order unity relating the observed bandpass to the standard 

restframe-bandpass, and the opening angle 9 jet is a function of the jet break time tbrcak = ijct. -E'iso(7), and 
the ambient number density n (see equation 1 of FKS). Since 9 jet is implicitly a function of Di, z, k, and 
S'-y, the value for E^ is a complex function of observables. 

Expressing all physical quantities in cgs units and making the analogy with astronomical magnitudes 
(distance modulus = 51ogio[-D;/10 pc]), we rearrange equation 1 to give, 
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where the DM is the "apparent GRB distance modulus." Assuming that the total energy output is constant 
(f-y) from burst to burst, equation 2 becomes. 
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with the zeropoint zp = 30.4940 + 2.5 logio(£7/l-5 foe). The unit "foe" equals 10^^ erg 



The variance of E^ is a direct measure of the accuracy with which E^ can be used for cosmographic 
purposes. The most direct method to estimate the variance on E~^ is to undertake a simultaneous analysis 
of the afterglow data and prompt burst emission to determine ft and E^^. With this direct approach we do 
not have to worry about hidden correlations between different quantities. 

Another approach is to carry out a simple error propagation of equation 1. The approach is straightfor- 
ward but does require the assumption that there is no covariance in the measurement of Sj and the inference 
of ^jet- We justify this assumption by noting that fh is obtained from afterglow measurements whereas 
derives from the prompt emission - two very different phenomena. There could well be global correlations 
between the prompt burst emission and afterglow emission (e.g., the brightest bursts may be associated with 
bright afterglow) but we see no reason why our inference of fi, should be observationally correlated to S^, 
beyond the small mathematical dependence of fb upon E^; see equation 1 of Frail et al. 2001). Adopting 
this approach toward variance estimation, we find the fractional error on the E^ measurement from. 
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where, 



9jet Sin^jet 



(5) 



_8(l-C0Sejet). 

The error in the ^-correction (cTfe) is given by equation 8 of Bloom et al. (2001b) and the error on the 
fluence (ag^) is estimated from the prompt GRB observations. The errors on density (fT„) and the jet- 
break time ((7/,j^, ) are foimd from afterglow modeling. Again, we have assumed that the measurements of 
the four observables are uncorrelated'^ . The error on the apparent GRB DM is 2.5 xas^/E^- If ^jet 1) 
then Cg.^^ « 1/16, implying that the fractional error on Ej can be easily dominated by uncertainties in the 
fc-correction and Sj. 



4. Results 

In FKS we adopted a value of n = 0.1 cm~"^ for the circumburst density for all bursts. One clear step 
toward improving the measurements is a more realistic estimate of n from broadband modeling of the 
afterglow light curves. We have now included the best density determinations that exist for about 1/3 of 

the GRBs in Table 1. In the past, large systematic differences in n have been derived owing to incomplete 
datasets and the use of approximate relations for estimating afterglow parameters (Granot & Sari 2002; Frail 
et al. 2003). However, with more precise photometric data and the increasing sophistication of afterglow 
modeling, the true diversity of number densities is becoming a well-established observation. Modeling yields 
estimates in the range 0.1 cm^'^ < n < 30 cm^'^ and there is now little support for extremes of either high 
or low density. Herein we make use of published estimates of n and adopt a (new) canonical value of 10 
cm~^ for all remaining events. 



4.1. Energetics and a Standard Energy 

In Table 2, we show the computed energetics values for E~^ and the associated apparent GRB distance 
moduli, DM, for a cosmology with = 0.7, = 0.3, and Hq = 65 /igs kms~^Mpc~^. Of the 9 GRBs 
with measured z, tjet and n, the logarithmically-weighted mean total energy release is 1.16 foe ± 0.07 dex 
with a median energy of 1.33 foe. We therefore adopt a new standard energy of £-y = 1.33 h^^ foe with an 
error of ± 0.07 dex. This standard value is dependent upon the choice of cosmology. 

If we expand the sample to include all bursts with knowns redshifts and tjd but not density, the 
logarithmically-weighted mean energy is 1.21 foe ± 0.08 dex. That these two mean values agree suggests 
that the assignment of n = 10 cm~^ to the 7 bursts with only tjet measurements reasonably reflects the 
density of the population without n measurements. Figure 1 shows a histogram of E^ measurements, tightly 
clustered about S-^,. Note that the adopted value for now found using a proper inclusion of density into 
the E^ formulation, is a factor of 2.7 higher than that found in FKS. This is close to the analytic difference 
expected (« 3.2 for 9jet *C 1) by assuming n = 10 cm~^ instead of n = 0.1 cm~^, as in FKS. 

As seen in Figure 1, 11 of the 16 GRBs with E^ measurements fall within 5 a of £.y and five bursts 



^Thc determinations of tjot and n arc inferred from afterglow modeling and so are correlated to some degree. However, 
from the best observed bursts often tjet is found directly from optical/IR data whereas density is best constrained with radio 
measurements. Thus, we expect the correlation on tjet and n to be weak. Bloom et al. (2001b) discuss the possible correlations 
between 5^ and k. 
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(970508, 990123, 990510, 000418, 011211) are outside 5 a oi£^. The variance of all 16 GRBs with measured 
about our adopted standard energy is 0.35 dex and is dominated by the five outliers. Though there is no 
a priori reason to exclude these outlier bursts, the observed scatter of the more restricted sample without 
the outliers about is 0.13 dex, about a factor of 35%. 

4.2. £^ in a Cosmological Context 

The small apparent scatter of measured E-^ values is continued evidence for a standard energy release 
in prompt 7-ray emission, first noted in FKS. As we pointed out in Garching 2001, this observation suggests 
that GRBs may be used as calibrated standard candles, useful, in principle, for measuring the geometric 
parameters of the universe ("cosmography"). 

A Hubble diagram, apparent distance modulus versus redshift, is one way to visualize the standard 
candle relation, now where the distance in the vertical axis corresponds to magnitude (i.e., 2.5 mag is a 
difference of a factor of 10). In Figure 2, we present the apparent GRB distance moduli for the 24 GRBs 
where such measurements (or constraints thereon) are now possible. The solid curves show the theoretical 
distance moduli in six cosmologies. For each cosmology, we re-compute £^ using the median of the same 
well-studied bursts as described above. This value of E-y is then used to set the zeropoint of the apparent 
GRB DM measurements from equation 3. 

Aside from an overall zeropoint, the values of the apparent GRB DM are almost"^ entirely dependent on 
observables, and not cosmology, yet the theoretical DM curve are sensitive to cosmology (oc 51ogrf;). That 
the theoretical DM curves trace the apparent GRB DM measurements in Figure 2, is simply a recasting of 
the standard candle energy result. 

4.3. Limitations in Cosmographic Applications 

One striking feature of Figure 2 is that the data appear to be well-represented by the theoretical curves 

in all six (rather different) cosmologies, with no apparent discrepant trend with redshift. The rms scatter 
about the respective is about the same in each cosmology (« 0.35 dex). dominated by the outliers in 
the distribution. The main difference between various cosmological formulations is that the standard energy 
varies by a factor of ~ 3, all about energies of -^1 foe. 

We explain the insensitivity of the energy data to cosmographic discrimination as follows. Without an 
a priori knowledge of the true standard energy release either via observations of a local sample or from 
theoretical studies — the data themselves are used to find the standard energy, setting the zeropoint for 
the DM measurements. This is akin to side-stepping the cosmological distance "ladder." Indeed without an 
anchor, the divergence of theoretical DMs in different cosmologies at large redshifts is subtle, and hence more 
difficult to observe. In other words, unlike with other cosmological tools such as with SNe Type la (with 
local calibrations) or the Sunyaev-Zeldovich effect (which depend only on local physics), any cosmographic 



^There is a small dependence of 9jet upon cosmology, such that for small opening angles, the value of DM goes as 1.25 x 
log rD;(x), for a ratio, rDi(z), of luminosity distances versus redshift relation in two different cosmologies. Even at the 
extrema, for most plausible cosmologies, this ratio is less than about 2 for redshifts less than five. Therefore, at most, the DM 
constructed from observable parameters could be systematically off by < 0..35; that is, the role of cosmology in determining the 
RHS of equation 4 is usually small compared with the error on the observables themselves. 
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discrimination with GRBs is performed differentially. As can be seen in Figure 2, the current scatter in 
GRB energies is simply to large for any apparent differential trend with redshift to be seen. The detection of 
lower-redshift {z ^ 0.3) bursts with Swift would help pin down £^ and increase the sensitivity of the sample 
to cosmology. 

5. Toward a More Homogeneous Set 

As noted in §1, cosmographic applications benefit enormously if the objects are well understood phys- 
ically. Lacking a clear understanding (as is arguably still the case with Type la SNe), an astronomical 
appreciation of the diversity of the phenomenon is critical. In short, if the scatter on is ever to be 
reduced a priori, it is essential that wc develop the ability to identify sub-classes from observations so that 
a homogeneous set of GRBs can be identified for cosmographic applications. 

We already know of at least one member of a distinct sub-class of GRBs: the GRB associated with 
SN1998bw (Galama et al. 1998; Kulkarni et al. 1998). Indeed GRB 980425, with an isotropic-equivalent 
energy release of ~ 7 x 10^^ erg, is clearly abnormal (see Figure 3). Bloom et al. (1998) suggested 
that such low- luminosity GRBs that are associated with supernovac to be "S-GRBs". Recently, Norris 
(2002) identified ^^90 such low-redshift S-GRB candidates in the BATSE sample. Without appealing to 
7-ray properties, such GRBs are readily identified by their low redshift, and possibly, the dominance of the 
associated supernova emission over the afterglow light. 

Of the more classical "cosmological GRBs," those with known redshifts beyond ^ 100 Mpc or associated 
with distant galaxies, we can use Figure 3 help identify classes of outliers. In Figure 3, we show the E^ 
distribution energy versus redshift for the 24 bursts with known redshift. We also show the trajectory of 
Ej for the 5 bursts where no redshift is known. These curves were found by computing the individual 
fc-corrections and /(, as a function of redshift. Note that the implied Ej curve begins to flatten at high 
redshift as a result of the mitigation between larger £"180(7) and smaller inferred /t. From Figure 3, it can be 
seen that the energy release in GRB 980329 and GRB 000630 could have been consistent with £j for redshift 
2: ^ 1.5. The present observations provide only a weak lower-limit to the energy of GRB 981226. 

5.1. "fast-fading GRBs" (f-GRBs)— A New Sub-Classification? 

As is readily seen in Figure 3, GRB 980519 and GRB 980326, stand out for being significantly lower 
in energy release than the cosmological set. Even at a redshift of 5, these bursts released less than 0.2 
foe. Interestingly, independent of the inferred low energy, GRB 980519 and GRB 980326 have previously 
been recognized as peculiar in that both afterglows exhibited an unusually steep decline, flux oc with 
a ^ 2 (instead of the usual a ^ 1) at early times. We propose that bursts with early-time steeply declining 
afterglows may be intrinsically under-energetic and suggest a nomenclature, "fast-fading GRBs" or f-GRBs 
(in analogy with the "Branch peculiar" of Type la SNe), for such GRBs. 

While the identification of f-GRBs, based upon fast-declining afterglows, is also an appeal to the physical 
state of those systems, under the assumption of an ISM density, we cannot identify any unifying observation 

of the 5 closer energy outliers (970508, 990123, 990510, 000418, 011211) that belie some unifying physical 
process. We therefore cannot (yet) distinguish such bursts a priori from the remainder of the set. 

Bursts that are f-GRBs could be either genuinely low-energy events or GRBs which typical energy 
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releases by taking place in particularly dense media (e.g., Dai & Lu 2001); note that for the best studied 
f-GRB — 030329 (see below) — there is no indication of a dense medium. In formulating our energy measure- 
ments we have explicitly assumed that the afterglow emission before the jet break arises in a medium where 
density is constant and homogeneous (i.e., p{r) oc r°). Since the environments of a least some GRBs are 
thought to be clumpy (Lazzati et al. 2002; Bcrsier et al. 2003; Heyl & Pcrna 2003) and wind-blown (i.e., 
p{r) oc r~^; Chevalier & Li 1999; Price ct al. 2002b), this is clearly not the case for all GRBs in our sample; 
therefore, our estimations of may be inappropriate for some GRBs. As in the homogeneous case, GRBs 
the dynamics of F in jetted GRB within winds should be insensitive to the jet edge so long as F <^ 1/^jet- 
We can therefore calculate ^jet(wind) by inverting the equation for F in eq. 1 of Chevalier & Li (2000): 

,.,wi„.,.«.0.(l±i)-"*(|Sa)-'\,..(i|)-'\ 

For all bursts in our sample for which A* has been reported we calculate 6'jet(wind) following cqs. 1 and 
6. The results are as follows (6'jet(wind), E^/IQ^° erg. A* reference): 970508 (20.81°, 3.560; Panaitescu & 
Kumar 2002); 980519 (4.88°, 6.262; Panaitescu & Kumar 2002) (5.13°, 6.941; Chevalier & Li 1999); 991208 
(<5.29°, <4.771; Li & Chevalier 2001); 990123 (2.82°, 17.390; Panaitescu & Kumar 2002); 000301C (8.15°, 
4.411; Li & Chevalier 2001); 000418 (11.78°, 15.834; Panaitescu & Kumar 2002); 000926 (5.24°, 11.689; 
Panaitescu & Kumar 2002); 010222 (1.93°, 4.891; Panaitescu & Kumar 2002); 011121 (>4.47°, >1.390; 
Price et al. 2002b). (We assume z = 2.5 for GRB 980519. Note that the wind hypothesis for GRB 991208 is 
strongly disfavored with late-time radio modeling of Galama et al. 2003a.) Comparison with Table 2 shows 
that the inferred values for and E-^ arc usually smaller under the wind assumption. Uncertainties in A* 
are not generally provided so it is unclear if the energy determinations are consistent between winds and ISM 
models. GRB 970508 is underluminous by a factor of ~3 in both the ISM and wind cases. Both GRB 000418 
and GRB 990123, overluminous in the ISM case, are more consistent with the standard £-y in the wind case. 

Only GRB 980519 and GRB 011121 show convincing evidence for a wind-stratified medium under the 
assumption of a top-hat jet model (GRB 980326 is consistent with a wind; Bloom ct al. 1999). As seen, the 
energy release in GRB 011121 is similar in both the wind and ISM cases. Interestingly, the energy release 
in GRB 980519, again assuming z = 2.5 and a wind-blown circumburst medium, is 6.3 xl0^° erg and thus 
consistent with f^. Therefore it may be reasonable to suspect that all f-GRBs are bursts that occur in 
wind-stratified media (but not necessarily the converse). 



6. Discussion and Conclusions 

With nine new bursts suitable for measurement after FKS, the evidence for a standard release of 
energy in GRBs persists. By incorporating more realistic measurements of ambient density and under the 
assiimption of a top-hat jet and constant ambient density media, wc find a standard candle energy release 
of £^ = 1.33 /igg^ foe, with a burst-to-burst scatter of 0.35 dex about this value. Though this result makes it 
interesting to consider cosmographic applications of GRBs, we have shown that without a local calibration 
of the true energy release, the observed scatter is simply too large to make any meaningful impact on 
cosmographic measurements. 

The most obvious way to reduce the scatter on £^ is to simply remove outliers from the sample, based 
purely upon discrepant energy measurements. However, this is a dangerously circular approach since it is 
particularly these apparent outliers in energy, at high redshift especially, that place the most meaningful 
constraints for cosmography. 
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Still, at this stage, we are confined to first recognize heterogeneity in the context of a cosmological 
model and then find similarities between discrepant bursts, based upon observations or theory that are 
complementary to energy. To this end, we have noted that both GRB 980519 and GRB 980326 are low- 
energy bursts and stand out based upon afterglow observations. 

Indeed, lacking firm theoretical predictions for sub-classification, we have proposed that there are at least 
two sub-classes - S-GRBs and f-GRBs - of GRBs which can be typed or classified from simple observations. 
Both groups were first identified by low inferred 7-ray energy releases. Now, S-GRBs arc characterized by 
low redshift^ and f-GRBs are identified by the rapid fading, of their afterglow at early times, t < 0.5 
day (restframe). 

The advantage of empirical distance indicators (e.g., Schaefer 2003) is that they can be entirely based 
on the 7-ray observations and do not require (once calibrated) additional time-consuming observations from 
the groimd. However, as with many empirical indicators, the method relics on poorly understood physical 
mechanisms (although, see, Salmonson 2000). The sheer bewildering diversity in any high-energy property 
of GRBs - richness of profiles, overall profile shapes, fluences, peak luminosities and even peak energies - 
should give pause to claims of the predictive power of distance indicators in the absence of a realistic physical 
model for the emission mechanisms. 

The physical uncertainty aside, it is clear that the high energy properties can and have provided gross 
distinctions. The well-known separation in energy and spectra between short-bursts and hard-bursts (Kouve- 
liotou et al. 1993), for example, has led to the belief that bursts from each class may originate from different 
progenitors. As a more recent example, we note that Norris (2002) has used the lag Imninosity relation to 
identify a subsample of long-duration bursts whose properties distinguish them from the general popula- 
tion. However, we question whether empirical relations have demonstrated that they possess the precision 
necessary to advance cosmography. There is little reason based on the limited understanding of the GRB 
phenomenon to strongly justify an optimistic expectation in this direction. 

The approach we have used in FKS and herein has the advantage that it is based upon energy release, 
a fundamental quantity in explosions; on physical grounds, it is reasonable to expect that such a quantity 
could be standard (or at least amenable to calibration) despite the rich diversity in GRB properties. The 
standard energy result is corroborated by other studies (which do not rely on the 7-ray data), namely multi- 
wavelength afterglow modeling and X-ray observations (Panaitescu & Kumar 2001; Piran et al. 2001; Berger 
et al. 2003). 

We stress that if GRBs are ever to be used for cosmographic purposes, there must be a significantly 
increased understanding in the diversity of the phenomenon (let alone the secular evolution with redshift). 
To this end, using simple observations, at least two sub-classes of GRBs, the S-GRBs and the f-GRBs, may 

be readily recognized. For the short term, we believe that the identification of such bursts will provide infor- 
mative direction to hone in on GRB progenitors. To weed out such bursts and create a more homogeneous 
sample for cosmography, however, our proposed identification methods have the disadvantage of requiring 
significant follow up observations. Nevertheless, it is possible that in the SWIFT era., dedicated facilities will 
undertake the necessary follow-up observations quite routinely. 

After this paper was submitted, GRB 030329 was discovered and found to fade rapidly at early times 



''The origin of 7-ray emission in S-GRBs has now been discussed in some detail (e.g., Tan et al. 2001) and it is reasonable to 
infer that the 7-ray emission arises from entirely different process as compared to the internal shocks of cosmologically located 
GRBs. 
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{F^ oc i;i-9'i'±o.i2-j (■pj.jgg oi 2003a) suggesting that the burst fit our definition of an f-GRB. Using the same 
formahsm above, we find = (8.18 ± 1.00) x 10**^ erg assuming a constant density ISM and n = 10 cm~^ 
or 1.3 X 10^° erg assuming a wind with A* = 1.0. In the former case, GRB 030329 is underluminous by a 
factor of 16.3 from the £^ found herein. Thus GRB 030329 joins GRB 980519 and GRB 980326 in the same 
sub-class both on energetics and phenomenological grounds. 

Given new the connection of GRB 030329 with a supernova (Stanek ct al. 2003), we suggest that all 
f-GRBs may be associated with supernovae. Obviously with afterglows that fade fast there is an increasing 
likelihood of detecting a supernova component (the non-detection of a SN component in the optical light 
curve of GRB 980519, see Jaunsen et al. 2001a, may simply be due to 2; ^ 1 of this burst). However, we 
suggest that the connection may be more than due to observational bias; instead long-duration GRBs could 
arise from two populations of bursts (Chevalier & Li 1999), with f-GRBs from massive stars that produce 
winds prior to the GRB and a concurrent supernova with the GRB event. Alternatively, especially in light 
of the detection of a supernova component with GRB 011121 (Bloom et al. 2002; Garnavich et al. 2003), we 
favor the hypothesis that all long-duration GRBs also lead to supernovae, with the diversity in GRB energies 
set by a progenitor-dependent ratio of the energy channeled into the relativistic jets to the energy powering 
the explosion the massive star progenitor. 
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Fig. 1. — A histogram of GRB energies (E^) with three equal logarithmic spacings per decade. The histogram 
shows a narrow distribution of GRB energies about the standard energy £^ = 1.33 foe, with an error of 
(J = 0.07 dex. The observed burst-to-burst rms spread is 0.35 dex (a factor of 2.23) about this value. Bands 
of 1, 2, and 5 a about the standard energy are shown. There are five identifiable outliers, which lie more 
than 5 a from the mean, however, there is currently no basis other than discrepant energy to exclude these 
bursts from the sample. In Figure 3 we identify two bursts (not shown here) which are discrepant in both 
energy and afterglow properties. 
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The GRB Hubble Diagram with 24 sources 
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n^ = 0.30, ^^=0.70 
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Fig. 2. The GRB Hubble Diagram illustrating the difficulty in making cosmographic measurements with GRBs. Shown 
are Hubble Diagrams for six different cosmologies for the 24 GRBs with known redshift and a measurement or constraint on 
tjet- The solid squares represent the apparent GRB distance modulus in the given cosmology with associated 1 a errors. Open 
circles show those sources with upper or lower limits on the DM measurement. The solid curves are the theoretical DM for the 
given cosmology. Though the values of Ey vary by more than a factor of 3 over the cosmologies shown, the observed variance 
about the theoretical curves is almost the same (Ri 0.35 dex); we discuss in the text why this is so. Clearly, the current sample 
cannot distinguish between a wide range of cosmologies. Here, and throughout the paper, we use Hq = 65 km Mpc~^. 
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0.30 0.4 0.6 1 2 3 4 5 6 

redshift (z) 

Fig. 3. — GRB energy release versus redshift assuming a top-hat jet illustrating the diversity in GRB energies. 
Bands of 1, 2, and 5 a about the mean energy = 1.33 foe are shown. Solid box and open circle symbols 
are the same as in Figure 2. Plotted are the trajectories of five GRBs with no known spectroscopic redshift 
(labeled with star symbols). Redshift limits are place from photometric observations of the afterglow (980519: 
Jaunsen et al. 2001a; 980326: Bloom et al. 1999; 000630: Fynbo et al. 2001; 991216: Halpern et al. 2000) 
assuming no extinction from the redshifted Lyman a forest in the most blue filter. The photometric redshift 
of the host galaxy of 980329 is from Jaunsen et al. (2003). While the energies of GRB 980329 and GRB 000630 
could be consistent with the standard value at redshifts beyond z ~ 1.5, at no redshift could the energies 
of GRB 980326 and GRB 980519 (dashed lines) be consistent unless the densities were significantly higher 
than the canonical value of 10 cm^''. Interestingly, these two bursts are associated with having occurred in 
a wind-blown density environment. While most long-duration GRBs appear to fall within a narrow range of 
energies, this shows there are several exceptional outliers, with demonstrably different afterglow properties 
from other bursts; we refer to such bursts as "fast-fading GRBs" or f-GRBs. 
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Note. — References are given in order for the redshift, fiuenee, jet-break time (ijct): and density (n). Where several redshift 
measurements are available for the same burst, we choose the most precise determination. For fluence measurement, we choose the 
measurement with the most precisely measured prompt burst spectrum. For bursts with more than one tjet or n measurement, we 
choose the reference where the most broadband data were used for afterglow modeling, preferring modeling that includes radio afterglow 
measurements. 

*The energy range over which the fluence was reported. 

^The prompt burst spectral /c-correction, as defined in Bloom et al. (2001b), used to transform the observed fluence in the 
particular bandpass to a standard restframe bandpass of 20—2000 keV. 

References. l. Bloom ct al, (2001a); 2. Frontcra ct al. (1998); 3- Bloom ot al. (199S); 4. Piran ct al. (2000); 5. Frail ct al. (2nOnb): 6. Djorgov.ski 

et al. (2001); 7. Kulkarni et al. (1998); 8. Groot ct al. (1998); 9. Yost ct al. (2002): 10. Galama ct al. (1999): 11. Jaunscn ct al. (2001b); 12. Paiiaitcscu & 
Kumar (2002); 13. Djorgovski et al. (2000); 14. Woods ct al. (1998); 15. Halpcrn Sz Fcscn (1998); 16. Djorgovski ct al. (1998); 17. Frail ct al. (2003); 18. 
Prontera et al. (2000); 19. Frail et al. (1999); 20. Kulkarni et al. (1999); 21. Bloom et al. (2003); 22. Vreeswijk et al. (2001); 23. Harrison et al. (1999); see 
also Stanek et al. (1999); 24. Le Floc'h et al. (2002); 25. Amati et al. (2000); 26. Masetti et al. (2000); 27. Frontera et al. (2000); 28. Fruchter et al. (2000); 
29. Castro-Tirado et al. (2001); 30. Hurley et al. (2000); 31. Sagar et al. (2000); 32. Vreeswijk et al. (1999); 33. Halpern et al. (2000); 34. Andersen et al. 
(2000); 35. Piro et al. (2002); 36. Castro et al. (2000b); 37. Berger et al. (2000); 38. Berger et al. (2001); 39. Pynbo et al. (2001); 40. Price et al. (2002a); 
41. Castro et al. (2000a); 42. Price et al. (2001); 43. Harrison et al. (2001); 44. Mirabal et al. (2002); 45. in' t Zand et al. (2001); 46. Galama et al. (2003b); 
see also Jha et al. (2001); 47. Price et al. (2002c); 48. Ricker et al. (2002); 49. Price et al. (2002); 50. Garnavich et al. (2003); 51. Price et al. (2002b); 52. 
Holland ct al. (2002); 53. Frontcra ct al. (2001); 54. Fox et al. (2002); 55. Ricker et al. (2002a); 56. Price et al. (2003b); 57. Ricker et al. (2002b); 58. Barth 
ct al. (2003); 59. Hurley ct al. (2002); 60. Matheson et al. (2002); 61. Lamb et al. (2002); 62. Pandey et al. (2002); 63. Schaefer et al. (2002); 64. Vreeswijk 
et al. (2002); 65. Crew et al. (2002) 
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Table 2. Computed GRB Energetics and the Apparent GRB Distance Modulus 
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> 9.84 ± 1.47 


< 47.83 ± 0.37 


0.26 


0.18 


0, 


,06 


0, 


,19 


980613 


53.55 ± 10.39 


> 12.57 ± 0.58 


> 1.28 ± 0.28 


< 47.03 ± 0.54 


0..38 


0.33 


0, 


,06 


0, 


,19 


980703 


601.20 ± 65.54 


11.21 ± 0.81 


11.47 ± 2.07 


44.31 ± 0.45 


0.26 


0.11 


0, 


,22 


0, 


,27 


990123 


14378.92 ± 1778.12 


4.93 ± 0.43 


53.14 ± 11.33 


44.00 ± 0.53 


0.26 


0.19 


0, 


,06 


0, 


,42 


990506 


8614.96 ± 1012.56 




















990510 


1763.49 ± 200.18 


3.36 ±0.21 


3.04 ±0.51 


47.13 ±0.42 


0.26 


0.14 


0, 


,28 


0, 


,12 


990705 


2559.52 ± 203.18 


5.33 ±0.41 


11.05 ± 1.91 


43.98 ± 0.43 


0.06 


0.20 


0, 


,06 


0, 


,37 


990712 


49.70 ± 5.27 


> 40.80 ± 1.70 


> 12.08 ± 1.60 


< 42.13 ± 0.33 


0.12 


0.25 


0, 


,06 


0, 


18 


991208 


1121.55 ± 329.97 


< 8.64 ±0.77 


< 12.72 ±4.38 


> 43.35 ± 0.86 


0.26 


0.71 


0, 


,36 


0, 


19 


991216 


5353.69 ± 593.92 


4.57 ±0.72 


17.03 ± 5.65 


44.02 ± 0.83 


0.26 


0.12 


0, 


,46 


0, 


,62 


000131 


11618.99 ± 1767.35 


< 4.68 ±0.21 


< 38.66 ± 6.79 


> 47.06 ± 0.44 


0.26 


0.30 


0, 


,06 


0, 


,19 


000210 


1693.26 ± 140.51 


> 6.84 ±0.28 


> 12.04 ± 1.40 


< 43.90 ± 0.29 


0.08 


0.20 


0, 


,06 


0, 


19 


00030 IC 


437.49 ± 62.21 


13.14 ±0.51 


11.46 ± 1.86 


46.30 ± 0.40 


0.26 


0.26 


0, 


,12 


0, 


13 


000418 


751.37 ±76.31 


22.30 ±6.34 


56.20 ±32.06 


42.97 ±1.43 


0.26 


0.05 


1, 


,35 


0, 


,37 


000911 


3954.79 ±630.19 


< 5.63 ± 0.25 


< 19.09 ± 3.48 


> 44.00 ± 0.46 


0.26 


0.32 


0, 


,06 


0, 


,19 


000926 


2797.36 ± 990.19 


6.16 ±0.31 


16.17 ±5.95 


45.93 ± 0.92 


0.26 


0.87 


0. 


,07 


0. 


,10 


010222 


8578.41 ± 216.66 


3.20 ±0.13 


13.33 ± 1.13 


45.28 ± 0.21 


0.06 


0.01 


0, 


,06 


0. 


,19 


010921 


136.11 ± 26.70 


31.19 ±2.46 


19.67 ±4.90 


41.71 ±0.62 


0.03 


0.50 


0, 


,06 


0, 


,36 


011121 


455.59 ± 188.54 


> 15.35 ± 1.00 


> 16.26 ± 7.05 


< 41.35 ± 1.08 


0.26 


1.03 


0, 


,06 


0, 


19 


011211 


672.34 ± 85.68 


6.38 ±0.40 


4.17 ±0.74 


47.54 ±0.45 


0.26 


0.20 


0, 


,06 


0, 


,30 


020405 


720.08 ± 92.23 


7.81 ±0.93 


6.68 ± 1.80 


43.99 ±0.67 


0.27 


0.19 


0, 


,06 


0, 


,58 


020813 


7749.76 ± 3105.50 


3.13 ±0.23 


11.58 ±4.94 


45.00 ±1.07 


0.26 


1.00 


0, 


,06 


0, 


,26 


021004 


556.01 ±72.10 


12.67 ±4.51 


13.53 ± 9.77 


46.48 ± 1.80 


0.26 


0.21 


1, 


,77 


0, 


,07 


021211 


378.24 ± 201.58 





















^E^ is computed using equation 1. The errors on the four observable parameters k, tjet, S-,, and n are taken from Table 1. 
When no errors are known for a given parameter we assume a fractional error of 10%. We have taken = 0.3 and S^a = 0.7 
for these energetics calculations. 

''DM is apparent GRB distance modulus calculated from equation 3, assuiniug f,, = 1.33 foe. The cosmology as noted 
above was used for the calculating, but note that the measurement on DM is relatively insensitive to the choice of cosmological 
parameters (see footnote 3). 

'^These columns give the error contribution to the DM measurement from each observable parameter. The error on DM is 
the quadrature sum of these four numbers. Divide these numbers by 2.5 to get fractional error contribution to E-^. 
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